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ABSTRACT 

This paper is the third in a series in which we present deep Chandra ACIS- 
S imaging spectroscopy of the Seyfert 1 galaxy NGC 4151, devoted to study 
its complex circum-nuclear X-ray emission. Emission features in the soft X-ray 
spectrum of the bright extended emission (Xo.3-2keV ~ 10 40 erg s _1 ) at r > 
130 pc (2") are consistent with blended brighter OVII, OVIII, and NelX lines 
seen in the Chandra HETGS and XMM-Newton RGS spectra below 2 keV. We 
construct emission line images of these features and find good morphological 
correlations with the narrow line region clouds mapped in [OIII]A5007A. Self- 
consistent photoionization models provide good descriptions of the spectra of the 
large scale emission, as well as resolved structures, supporting the dominant role 
of nuclear photoionization, although displacement of optical and X-ray features 
implies a more complex medium. Collisionally ionized emission is estimated to 
be <12% of the extended emission. Presence of both low and high ionization 
spectral components and extended emission in the X-ray image perpendicular to 
the bicone indicates leakage of nuclear ionization, likely filtered through warm 
absorbers, instead of being blocked by a continuous obscuring torus. The ratios 
of [OIII]/soft X-ray flux are approximately constant (~15) for the 1.5 kpc radius 
spanned by these measurements, indicating similar relative contributions from 
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the low- and high-ionization gas phases at different radial distances from the 
nucleus. If the [OIII] and X-ray emission arise from a single photoionized medium, 
this further implies an outflow with a wind-like density profile. Using spatially 
resolved X-ray features, we estimate that the mass outflow rate in NGC 4151 is 
~2M yr _1 at 130 pc and the kinematic power of the ionized outflow is 1.7 x 10 41 
erg s" 1 , approximately 0.3% of the bolometric luminosity of the active nucleus 
in NGC 4151. 

Subject headings: X-rays: galaxies — galaxies: Seyfert — galaxies: jets — galax- 
ies: individual (NGC 4151) 

1. Introduction 

The Seyfert 1 galaxy NGC 4151 is well known to have an extended narrow line region 
(ENLR) in the optical (e.g., Perez et al. 1989; Robinson et al. 1994; Crenshaw et al. 2000), 
indicating interaction between the active galactic nucleus (AGN) and the host galaxy (Unger 
et al. 1987). Thanks to its proximity (d ~ 13.3 Mpc, 1" = 65 pc; Mundell et al. 1999) and 
extensive coverage at all wavelengths (Ulrich 2000 and references therein), NGC 4151 thus 
represents an excellent laboratory to study observationally nuclear feedback. 

Presence of extended emission in the soft X-rays is also well established (Elvis et al. 1983; 
Morse et al. 1995; Ogle et al. 2000; Yang et al. 2001; Wang et al. 2010a,2011a). The nature 
of the kpc-scale extended X-ray emission remains uncertain; suggestions include collisionally 
ionized plasma, photoionized gas, a hybrid of both, or electron-scattered nuclear emission 
(Elvis et al. 1983, 1990; Heckman & Balick 1983; Morse et al. 1995; Schulz & Komossa 
1993; Ogle et al. 2000; Komossa 2001; Yang et al. 2001), although electron-scattering is 
less favored, because the required scattering column and spectral variability do not match 
the observations (Weaver et al. 1994; Yang et al. 2001). Previous Chandra imaging studies 
of NGC 4151 (Ogle et al. 2000; Yang et al. 2001; Gonzalez Martin 2008; Wang et al. 
2009, 2010a) found good morphological agreements between the extended X-ray emission 
and the optical forbidden line emission (in particular [OIII] A5007). Such a close spatial 
correspondence between the soft X-rays and the [OIII] emission appears to be common in 
nearby Seyfert 2 galaxies, leading Bianchi et al. (2006) to propose that optical and X-ray 
features arise from a single photoionized medium. 

High spectral resolution grating observations of NGC 4151 with Chandra and XMM- 
Newton further show that the soft X-ray emission is due to both blended emission lines 
and radiative recombination continua (RRC) from He-like and H-like transitions of carbon, 
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oxygen, neon, and nitrogen (Ogle et al. 2000; Schurch et al. 2004; Armentrout et al. 2007). 
The X-ray line ratio diagnostics and electron temperatures of the RRC features support the 
conclusion that most of the gas is photoionized (T ~ 10 4 K) by the AGN (Schurch et al. 
2004; Armentrout et al. 2007), although photoexcitation by the AGN continuum may also 
play an important role (Kinkhabwalla et al. 2003, Ogle et al. 2003). 

However, these grating studies of NGC 4151 cannot provide the spatial distribution of 
each spectral feature. The extracted spectra consist of both nuclear emission and extended 
emission from various physical scales. Therefore, based on the high resolution spectra alone, 
one cannot exclude that the photoionized emission may be associated with the bright narrow- 
line region (NLR) gas clouds close to the nucleus (^C 100 pc), while a significant contribution 
from collisionally ionized gas may be present at larger radii (e.g., NGC 1365, Wang et al. 
2009, Guainazzi et al. 2009). Although Yang et al. (2001) presented ACIS imaging data, the 
limited signal to noise ratio (S/N) of these data prevented detailed high spatial resolution 
comparisons of optical [OIII] and X-ray morphology, or the spectral analysis of X-ray emission 
from the spatially resolved features. 

The X-ray emission mechanisms in this prototypical Seyfert 1 galaxy deserve further 
investigation with high quality X-ray data that allow spatially resolved spectral analysis of 
the circum-nuclear region. This is the aim of the study presented in this paper, the sixth 
and last in our series of studies on the nuclear and circum-nuclear region of NGC 4151 using 
deep Chandra ACIS and HRC observations (PI: Fabbiano). Using the HRC data and image 
deconvolution techniques, Wang et al. (2009) resolved the nuclear emission on spatial scale 
of ~30 pc and made comparisons with HS T observations. Wang et al. (2010a) presented the 
X-ray spectral analysis of the NGC 4151 nuclear emission. Wang et al. (2010b) reported the 
discovery of faint soft X-ray emission extending 2 kpc out from the active nucleus and filling 
in the cavity of the H I material. In Paper I of the present series (Wang et al. 2011a), based 
on in-depth analysis of the ACIS data set, we have examined the kpc scale X-ray morphology 
and the relations between X-ray absorption features and the cold ISM in the host galaxy. In 
Paper II (Wang et al. 2011b) we have reported on strong evidence of jet-cloud interaction 
in the nuclear region (r < 130 pc) and compared these findings to our previous HRC study. 
Here we focus on the comparison between the X-ray line emission of the ionization cone and 
the ionized gas traced by the optical line emission, and we report on the spectral analysis of 
spatially resolved X-ray features. The paper is organized as follows. In § 2, we first briefly 
describe the data reduction, present the X-ray emission line morphologies (§ 2.1) and the 
spectral modeling of all extended emission within a 30"-radius (~ 2 kpc) and the spatially 
resolved emission (§ 2.2). The results are discussed in § 3 and summarized in § 4. 
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2. Data Extraction and Analysis 

NGC 4151 was observed by Chandra for a total of 180 ks (after screening for high back- 
ground intervals) with the spectroscopic array of the Advanced CCD Imaging Spectrometer 
(ACIS-S; Garmire et al. 2003) in 1/8 sub-array mode during March 27-29, 2008. We have 
presented the details of our ACIS observations of NGC 4151 and the data reduction in Paper 
I. 

The data were analyzed following the standard procedures using CIAO (Version 4.2) 
with the CALDB 4.2.1 provided by the Chandra X-ray Center (CXC). Subpixel event repo- 
sitioning and subpixel binning techniques (Paper I and references therein) have been applied 
to the ACIS images to improve the spatial resolution. 

The complexity in the data analysis caused by the bright nuclear emission was described 
in Papers I and II. We established that for the soft X-ray emission (E < 1 keV), photon 
pile-up is most severe in the inner r < 1", mild at 1" < r < 2", and not an issue at 
r > 2". Chandra PSF simulations were performed to provide an estimate of the expected 
contamination from the nuclear emission in an extended feature. We have taken into account 
this information in the following analysis. 

2.1. Broad-band Images and X-ray Emission Line Maps 

Figured] (top) shows the 0.3-2 keV ACIS image (~0.5" per pixel). Bright structured soft 
X-ray emission along the northeast (NE) - southwest (SW) direction is apparent, together 
with fainter, less extended emission along the northwest (NW) - southeast (SE) direction. 
The inner r < 2" nuclear region, where we identified interactions between the radio jet 
and the optical clouds, has been discussed in Paper II. Here we focus on the kpc-scale 
extended emission. We first extracted counts from every 10° sectors to create azimuthal 
surface brightness profiles for the extended X-ray emission (2" < r < 30"), which provides a 
clear view of its angular distribution, and to identify sectors that contain the bright outflows 
(Figure [U bottom). This plot was used to guide further spectral extraction (see § I2.2.2p . 
The extended emission along the NW-SE direction is interesting as this is the direction of 
the putative obscuring torus. Figure [2] demonstrates the clear elongation perpendicular to 
the bicone (in particular the NW sector), using the 0.3-2 keV ACIS image smoothed with 
the CIAO tool csmooth. We extracted the 0.3-2 keV radial profiles of the NW sector, the 
bright SW cone, and the faint emission in between (the "control region"), shown in Figure [3J 
The NW emission is on average 5 times fainter than the bright SW cone, but brighter than 
the control region at r < 7" (4cr significance). Beyond r = 7", it becomes indistinguishable 
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with the faint emission. The implication on the nuclear obscuring structure will be discussed 
later in § IO 

A quick examination of the spectrum of the extended emission extracted from a circular 
region between a radius of 2" and 30" shows the clear presence of strong blended emission 
lines (Figure HJ), corresponding to the emission lines (<2 keV) seen in HETG observations 
(Ogle et al. 2000). These lines cannot be uniquely resolved with the spectral resolution of the 
ACIS imaging observation. Most notably the blended lines appear as three strong emission 
features approximately centered at 0.57 keV (OVII), 0.68 keV (OVIII), 0.91 keV (NelX) and 
a few weak lines between 1-2 keV (Mg XI, Si Ka, SiXIII; Figure Hj). The blended emission 
features seen in the ACIS spectrum and the corresponding HETG lines are summarized in 
Table 1 (see § 2.2.1). 

Following Paper II, we created adaptively smoothed images of the NGC 4151 circum- 
nuclear region using the subpixel resolution data. We extracted images from the merged 
data in three spectral bands dominated by emission lines below 2 keV: 0.3-0.7 keV ("soft 
band" containing the OVII and OVIII emission), 0.7-1.0 keV ("medium band" containing 
the NelX emission), and 1-2 keV ("hard band"). Figure Eh presents a false color composite 
image of the central ~ 45" x 45" (~3 kpc on a side), circum-nuclear region of NGC 4151, 
where the soft, medium, and hard band smoothed images are shown in red, green, and blue, 
respectively. Figure [5b zooms in to the central 1 kpc-radius region, emphasizing the bicone 
emission. For visualization of the faint features, the nuclear region (see Paper II) is saturated 
and excluded from this figure. The X-ray biconical morphology has been described in Paper 
I, and here we focus on the spectral differences. The extended bicone emission (> 2") appears 
rich in 0.3-1 keV emission (red and green). Similar color emission (mixed red/green) is also 
seen perpendicular to the bicone direction with a hint of harder (blue) emission. Clumps 
of 0.3-0.7 keV (red; OVII and OVIII) emission are also seen, some particularly prominent 
features are seen in the SW cone ~ 6" and ~ 10" from the nucleus. 

To further highlight the regions where the strong emission lines arise, we also extracted 
the X-ray emission in narrow energy intervals (0.53-0.63 keV; 0.63-0.73 keV; 0.85-0.95 keV) 
to create line strength images. This is a reasonable approach since the line emission domi- 
nates over the weak underlying continuum (Schurch et al. 2004) in these narrow bands. The 
resulting images are shown in Figure [61 The position of the nucleus is indicated with a black 
cross. OVII and OVIII line emission is prevalent in the extended regions. NelX emission, 
while extended along the general P.A. of the large scale extended emission, is bright in the 
central r = 2" but becomes fainter quickly at r > 3" compared to the OVII emission. We 
examined similar emission line maps for the weak lines in the 1-2 keV range (NeX, NelX 
RRC, MgXI) and found that they are highly concentrated within the inner r = 2" region. 
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This is better demonstrated in Figure [6H, where the radial profiles of the line emissions are 
shown. 



2.2. Spectral Analysis 



Spectra and instrument responses were generated using CIAO Version 4.2 and analyzed 
with XSPEC Version 12.6 (Arnaud 1996). Background spectra were taken from a source-free 
region from the same chip. Spectra were grouped to have at least 20 counts per energy bin 
to allow x 2 fitting. Unless otherwise noted, we restricted our modeling to photon energies 
above 0.3 keV, where the AC TS calibration is go ojj], and below 2 keV, above which the 
nuclear continuum dominates ( jWang et al.ll2010bl ). The 90% confidence interval for a single 
interesting parameter is reported for all fitting results. 



2.2.1. The 30" -radius Extended Emission as a Testbed 

We first characterize the X-ray spectrum of bright extended emission extracted from the 
2" < r < 30" circular region, exploring a phenomenological spectral fitting approach that was 
broadly adopted in previous studies at CCD resolution (Yang et al. 2001; see also Ogle et al. 
2000, Smith & Wilson 2001). The "continuum" in the ACIS spectrum (a combination of true 
underlying continuum and highly blended weak emission lines) was described by a smooth 
bremsstrahlung emission. We emphasize that the bremsstrahlung emission here should not 
be assigned any physical origin, whereas the derived line emission should be considered as a 
lower limit in such an approach. 

Repeating the technique used to fit the spectrum of the nucleus (Wang et al. 2010), 
narrow emission lines were then added to the model. The Chandra HETG and VMM-Newton 
RGS spectra (Ogle et al. 2000; Schurch et al. 2004; Armentrout et al. 2007) were used as 
a guide for the identification of the soft X-ray lines. The X-ray spectrum of the nucleus of 
NGC 4151 is notoriously complex and variable. The hard X-ray continuum is dominated by a 
power law component suppressed by a partially covering absorber and a Compton reflection 
component (Weaver et al. 1994; Schurch et al. 2004; De Rosa et al. 2007; Wang et al. 
2010b; Lubinski et al. 2010), an d the absorption con s ists of both a mildly-ion ized "warm 



absorber" and a cold absorber (IKraemer et al.l l2005t iGonzalez-Martm 1 120081 ). The PSF 



scattered contribution from the nuclear emission is taken into account in the spectral fitting, 



http : / / cxc . cf a . harvard . edu/cal/ Acis/ 



-7- 



using the nuclear co mponent fixed at the best fit nuclear spectral model and normalization 



( IWang et al.ll2010bf ) scaled to the expected fraction from PSF simulations. The extraction 
region contains 4.7% of the 0.3-2 keV nuclear emission due to PSF scattering, corresponding 
to ~25% of the extracted counts. This can be taken as the upper limit of contamination to 
the extended line emission, as the line emission dominates the nuclear spectrum. 

The identified emission line blends are summarized in Table [I] and the spectrum shown 
in Figure [7] (see also Figure H]). The absorption column required by the fit is consistent with 
the line-of-sight Galactic column towards NGC 4151 (N H = 2.1 x 10 20 cm -2 ; Murphy et 
al. 1993). The X-ray luminosity of the extended emission (between 130 pc and 2 kpc from 
the nucleus) corrected for the Galactic absorption is Lo.3-2keV = 1.1 ± 0.2 x 10 40 erg s _1 , 
compared with L 2 _iokeV = 2 x 10 42 erg s" 1 for the unabsorbed nuclear source (Wang et al. 
2010b). 

However, such a simplified model does not provide constraints on the emission mech- 
anism, and more physically meaningful spectral models are needed. Fitting the data with 
combinations of absorbed optically-thin thermal emission with solar abundances (APEC 
model; Smith et al. 2001) was attempted but gave poor results (reduced xl > 3). The abun- 
dances were then allowed to vary from solar abundance in the APEC model using the variable 
abundance APEC model (VAPEC), in which Zo, Z^ e , Z$i, and Zp e were left free while other 
elements were fixed at solar abundance. When Z is allowed to vary, the fit generally im- 
proves but still is far from satisfactory (xl ~ 2). Moreover, the required abundances become 
unphysically low (Z < O.O1Z ) as the fit attempts to reproduce a smooth continuum. 

To have a self-consistent model, we have made use of the Cloudy photoionization mod- 
eling code (Ferland et al. 1998). Although the resolution of the ACIS CCD spectrum 
prevents us from performing line-based diagnostics, e.g. deriving line ratios of the Helium- 
like triplets, spectral fitting with photoionization models still provides useful constraints 
about the photoionization status (as in Bianchi et al. 2006, 2010; Gonzalez-Martin et al. 
2010). Using Cloudy version C08.00, which enables a Cloudy/XSPEC interface (Porter et 
al. 2006), we modelled the < 2 keV X-ray spectrum assuming an open plane-parallel ge- 
ometry ("slab"). The dimensionless ionization parameter (Osterbrock & Ferland 2006) is 
defined as U = Q/(47rr 2 cn^), where i%h is the hydrogen number density, r is the distance to 
the inner face of a model slab, c is the speed of light, and Q = f^ 6eV L v jhv is the emitting 
rate of hydrogen ionizing photons (s _1 ) by the ionizing source. We assumed the broken 
power-law used by Kraemer et al. (2005) for the AGN continuum (see also Armentrout et 
al. 2007) . We varied U and N H , the column density of the slab to create spectral models 
covering a grid of parameters (—3 < logU < 3, 19.0 < log Nh < 23.5 cm -2 ), which were fed 
to XSPEC. 
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Armentrout et al. (2007) successfully modeled the XMM-Newton/RGS soft X-ray spec- 
trum of NGC 4151 with three photoionization components: a high ionization component 
(log [7 = 1.3, \ogN H = 23.0), a medium ionization component (log [7 = 0, logN H = 23.0), 
and a low ionization component (log [7 = —0.5, log Njj = 20.5). In our initial attempt, we 
included all these components, with the ionization parameters and the column densities fixed 
at these values and the normalizations set free. In addition, the photoionized emission was 
absorbed with a line-of-sight column density, N H>Lo s . The fit to the ACIS spectrum was poor 
(Xv — 6)- This is not a total surprise, given that these best fit parameters, are optimized for 
the brightest photoionized material close to the nucleus (10 _3 -1 pc; Armentrout et al. 2007) 
that dominates the soft X-ray spectrum of NGC 4151, whereas the emission of interest here 
arise on the r ~130 pc-1.5 kpc scale. 

When the (U, Njj) parameters for the photoionized components were set free to vary, 
we were able to obtain a statistically satisfactory fit (xl — 1-1) > reproducing the observed 
ACIS spectrum of the large scale extended X-ray emission well (Figured]). The line-of-sight 
column density required by the fit is N^in^ = 2 x 10 20 cm -2 , consistent with the Galactic 



column towards NGC 4151 ( jMurphy et al.lll996l ). Two components, a high ionization phase 



(log £7 = 0.8, logN H = 20.0) and a low ionization phase (log [7 = —0.25, log Nh = 19.4), are 
required to produce both the hydrogen-like and helium-like neon and oxygen line emissions. 
The Cloudy model predicts gas pressures of 6 x 10~ 7 dynes cm -2 and 1.3 x 10~ 7 dynes 
cm -2 for the high U component and low U component, respectively. If they are collocated, 
additional pressure-perhaps from a hot intercloud medium in a multiphase ISM (Elvis et al. 
1983, Ogle et al. 2000)-is needed to reach pressure equilibrium. Note that the log U = —0.25 
component closely resembles the "Low" model component in Armentrout et al. (2007) and 
the "D+Ea" X-ray absorption component in Kraemer et al. (2005). 

These model components still left small residuals at ~0.5 keV (the NVII Lya) and 
1.02 keV (the NeX Lya). Adding in a third photoionized component can further improve 
the fit (best fit parameters log U = 1.23, logA# = 19.2), but the improvement appears not 
significant when an F-test is performed (at 70% significance level). Thus the two-component 
model is used as a template for the following analysis of subregions. 

Given the satisfactory spectral fits with photoionized components and lack of clear 
signatures of collisional ionization in the ENLR, how much thermal emission from hot gas 
may be present in the circum-nuclear region of NGC 4151? We obtain a constraint on the 
possible presence of the hot ISM by invoking a collisionally ionized component in the best 
fit photoionization model. After adding a thermal emission component (APEC) to the 
photoionized emission, the fit is improved with a kT = 0.81 ±0.05 keV (at > 99% significant 
level by F-test) and the residual at 0.9-1.1 keV is mostly accounted for. This thermal 



-9 - 



component may contribute < 12% of the total soft X-ray emission, L t h,o.5-2keV = 6 x 10 38 
erg s _1 . As a cautionary note, although APEC component improves the fit statistically, it 
remains possible that the residuals are due to photoexcitation of Fe L-shell and the enhanced 
resonance lines (e.g., Ogle et al. 2003; Kinkhabwala et al. 2003) that are inadequately 
accounted for with the Cloudy photoionization model components. This will be explored 
with the high resolution HETG spectrum in future work. We will discuss the possible 
presence of the thermal component further. 

We emphasize that, because of likely complex geometrical dependencies and absorption 
between these components in the kpc scale emission, this modeling may still be over-simplistic 
to simultaneously reproduce the observed features in the broad band. 



2.2.2. Spectral Fitting of Selected Subregions 

To perform spectral analysis of the extended X-ray emission, the large scale region was 
divided into four sectors with position angles (P.A.) as outlined in Figure [H Extraction 
regions are termed the "SW", "NW", "SE" and "NE" regions respectively (Table 2; see 
Figured]). The bisectors of NE and SW correspond to the approximate P.A. of the ENLRs 
(e.g., Robinson et al. 1994). They are further divided into subregions along the radial 
direction (e.g., SW 1, 2, 3, and 4; NE 1, 2, and 3), with ~1000 counts (0.3-2 keV) in the 
fainter outer regions (SW4 and NE3). The NW and SE sectors cover regions (2" < r < 15") 
perpendicular to the bicone, which is the direction of the large scale weak bar (Asif et al. 
1998). 

Corresponding to the variations in color, spectral differences are clearly visible when the 
ACIS spectra across sectors are compared (Figure [9]). The fitting results are summarized in 
Table 2, and Figures [TO] [TT], and [121 show the corresponding spectra and the fits. We find 
that: 

1. Except for the outer- most regions (NE3 and SW4), a high ionization component is 
present in all regions (log [7 ~ 0.83-1.26). Lower ionization components are present 
along the bicone direction (NE-SW), with log [7 values ranging between -0.5 and 0. 
Recall that the fit to the integrated spectrum, which includes emission from all the 
subregions, finds the log £7 ~ 0.8 and log £7 ~ —0.25 components. 

2. In the NE3 and SW4 regions, a very low ionization component (log U > — 1) and a high 
column density (log TV// = 23.5) are required by the fit to improve the residuals in the 
Mgl Ka and Sil Ka. An F-test indicates that this fluorescent component significantly 
improve the fit at 98% confidence. 
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3. The second photoionized component in the NW-SE regions, perpendicular to the ion- 
ization cone, shows a lower ionization (logC/ ~ — 1). 

We also extracted spectra of the X-ray emission in the bicone regions that are prominent 
in 0.3-0.7 keV emission (or "OVII" rich, visually red in Figure [5]) and in 0.7-1 keV emission 
(or "NelX" rich, visually green in Figure |5]). The spectral fits are shown in Figure IT3"a.b 
respectively. Whereas the fitted ionization parameter for the "NelX" rich emission is close 
to the average value, the "OVII" rich emission is generally of lower ionization, likely due to 
higher number density. This is consistent with the appearance of clumpy material in the 
optical. 

3. Discussion 

3.1. Nature of ENLR Emission: Photoionization and Collisional Ionization 

Kiloparsec-scale extended X-ray emission is ubiquitously observed in nearby obscured 
Seyfert galaxies with ENLRs (Bianchi et al. 2006; Guainazzi et al. 2009; Bianchi et al. 
2010; Gonzalez-Martin et al. 2010), but the ionization mechanism for this emission has 
been a subject of active debate. For NGC 4151, Yang et al. (2001) proposed a two phase 
model for the extended X-ray emission, including a photoionized warm gas and a hotter 
collisionally ionized gas (see also Komossa 2001). The relatively low electron temperatures 

(T ~ 10 4 10 5 K) inferred from the radiative recombination continua (RRC; see Liedahl 

1999 for review) and the line ratios from the He-like triplets strongly support the presence 
of a photoionized component (Ogle et al. 2000; Schurch et al. 2004; Kraemer et al. 2005). 
However, does the observed soft spectrum entirely consist of photoionized emission, or does 
collisional ionization become more important than the photoionization mechanism beyond 
the nuclear region (Armentrout et al. 2007)? 

3.1.1. Evidences for Photoionization as the Dominant Ionization Mechanism 
Several of our results point to photoionization as the dominant mechanism. 

• Emission Line Gas Morphologies — First of all, the morphological agreement between 
soft X-ray and optical [OIII] emission is very strong. In particular, there is a good 
spatial correlation between [OIII] line emission and the three strong X-ray emission 
lines (Figures (6^-c). There is clearly an overall morphological coincidence with the 
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ionized gas bicone traced by [OIII] emission, as noted in previous X-ray studies (e.g., 
Ogle et al. 2000, Yang et al. 2001). A good agreement between the peaks of optical 
and X-ray line emissions in the inner r = 3" is found. 

• [OIII]-to-Soft X-ray Ratio of the HST -resolved Clouds — Bianchi et al. (2006) studied 
a sample of eight Seyfert 2 galaxies, and suggested that the same kpc-scale gas pho- 
toionized by the AGN continuum can simultaneously produce the observed X-ray and 
[OIII] emission. Our HRC and ACIS studies of the inner 2" (< 130 pc) emission in 
NGC 4151 (Wang et al. 2009 and Paper II) have shown that, for most of the clouds 
close to the nucleus, the ratios of [OIII] to soft X-ray emission are consistent with the 
typical ratios in Bianchi et al. (2006), the exception being clouds collisionally ionized 
by the radio jet. 

Our subpixel-processed ACIS images of NGC 4151 enable us to directly compare the 
extended X-ray emission and those of the resolved [OIII] clouds (Figure [6]). The main 
cloud features are labeled in Figure HH Using the calibrated continuum-subtracted 
HST/WFPC2 F502N image (Kaiser et al. 2000), we measured the [OIII] fluxes for the 
clouds and listed them in Table [3J ACIS spectra were extracted from the same regions 
and fitted with the continuum plus emission line model to obtain their 0.5-2 keV X-ray 
fluxes. Figure [15] shows the [OIII] to soft X-ray ratio for the distinct cloud features 
at various radii to the nucleus (~150 pc-1500 pc). The ratios listed in Table [3] are 
consistent with a constant [OIII]/X-ray flux ratio of ~15 for the 1.5 kpc range spanned 
by these measurements. The radially constant [OIII] /X-ray ratios indicate that the 
relative contributions of the low- and high-ionization gas phases at different radii are 
similar, as pointed out by our referee. If the [OIII] and X-ray emission arise from a 
single photoionized medium, as discussed in Bianchi et al. (2006), this further implies 
that the ionization parameter does not vary significantly to a large radial distance 
and the density is approximately decreasing expected for a freely expanding 

nuclear wind based on mass conservation (Kraemer et al. 2000; Bianchi et al. 2006). 

• Spectral Fitting — The spectral analysis shows that, a model with just two photoionized 
components of high- and low-ionization parameters can describe well the observed X- 
ray spectra in various regions across the circum-nuclear area, whereas in the spectral 
fits that invoke only collisionally ionized components strong emission line residuals 
persist. This supports the dominant role of nuclear photoionization. 

• Kinematics of the Ionized Gas — The full velocity field of the NLR for NGC 4151 
was mapped by Kaiser et al. (2000), which supports a hollow biconical outflow model 
for these clouds. At distances exceeding 4" (~260 pc) from the nucleus, Kaiser et al. 
(2000) find that the gas cloud kinematics are consistent with circular rotation in a 



- 12 - 



gravitational field, identical to the host ISM. Therefore, except for a few clouds in the 
nuclear region that appear disturbed by the impact of the radio jet (IKaiser et al.ll2000t 



Mundell et al.ll2003l see also Paper II), there is no evidence for strong kinematic shocks 



or any strong shocks in the outflow, which is consistent with the gas being photoionized 
by the nuclear radiation. 



Based on all the above considerations, we conclude that photoionization by the nucleus 
is most likely the dominant ionization mechanism in the ENLR of NGC 4151, confirming 
the suggestions of previous studies (Ogle et al. 2000; Schurch et al. 2004; Kraemer et al. 
2005; Armentrout et al. 2007; Wang et al. 2009b). It is reassuring to establish a role for 
nuclear photoionization in the extended NLRs that is consistent with the grating studies. 
However, there is also evidence that points to a more complex emission mechanism than a 
single photoionized medium (e.g., Bianchi et al. 2006). We discuss these results and their 
implications below in Sections 3.2 through 3.4. 



3.1.2. Constraints on the Presence of Collisionally Ionized Gas 

In § 12.21 we estimated the extent of thermal emission from hot gas (L < 6 x 10 38 erg 
s _1 ) that could be present in the circum-nuclear region of NGC 4151. The Chandra grating 
study of NGC 4151 by Ogle et al. (2000) provided the first direct evidence of X-ray line 
emission from a collisionally ionized hot plasma. Armentrout et al. (2007) also noted that 
the poor fit to the OVII line profile could be due to a non-photoionized component. Wang 
et al. (2011b) further show that a kT ~ 0.6 keV thermal emission component is present in 
the inner r = 2" region, which likely originates from jet-cloud interaction. 

Although star formation activity in the circum-nuclear region is at a low level (Asif et 
al. 2005; Dumas et al. 2010), hot gas may originate from the NGC 4151 galactic disk, since 
its temperature (kT ~ 0.8 keV) and luminosity (L ~ 6 x 10 38 ) are comparable to the disk 
diffuse emission in some star forming galaxies, kT ~ 0.6 keV, L ~ 10 39 erg s -1 (e.g., NGC 
253; Strickland et al. 2002). 

Another plausible alternative is that the hot plasma is the intercloud medium uniformly 
filling in the outflow (i.e., filling factor / = 1), providing pressure confinement for the NLR 
clouds (Elvis et al. 1983; cf. Krongold et al. 2007 for warm absorber outflows). The volume 
emission measure for the thermal component is n 2 e V = 4 x 10 61 cm -3 , derived from the 
normalization after fitting an additional thermal component (APEC) to the photoionization 
models of the bright biconical emission (§ 2.2.1). Assuming the hollow bi-conical NLR 
geometry (Crenshaw et al. 2000; Das et al. 2005), we obtain V ~ 10 61 cm -3 and so an average 
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density of ~ 2 cm -3 . Thus the thermal pressure of the T ~ 10 7 K gas is ~ 5.5 x 10~ 9 dyne 
cm -2 . Comparing with the thermal pressure of T ~ 10 4 K and n = 1600 cm -3 NLR clouds 
(Robinson et al. 1994) gives p = 4.8 x 10~ 9 dyne cm -2 . Thus pressure equilibrium between 
the hot and cool phase is highly likely, confirming the findings on pressure confinement in 
Ogle et al. (2000). 



3.2. A Complex Photoionized Medium: Displacement of Optical and X-ray 

Features 



Figure clearly shows that one distinct OVII emission blob "A" (at 6", ~390 pc SW 
from the nucleus) does not align with a bright [OIII] emitting cloud, designated as "Kl" in 
Robinson et al. (1994) and offset by 1" to the SE of "A" (see Yang et al. 2001). We find this 
is also the case for a more distant X-ray blob "B" (10" SW) where the closest [OIII] emission 
cloud is 2" to the SW. Similar displacements between the X-ray enhancement and the [OIII] 
blob are also clearly found in the OVIII line emission (Figure Eb). Note that the Chandra 
astrometry is accurate to 0.3"(Paper I), therefore the mis-alignment between features at this 
spatial scale is significant and indicates that the origin of the optical/X-ray emission needs 
further investigation. 

The radial location of X-ray blob "A" lies at the boundary between NLR and ENLR 
identified in the optical ( Robinson et al. 1994), where the density drop s rapidly from n ~ 1600 
cm -3 to n ~ 200 cm -3 (IPenston et al.l Il990l ; iRobinson et al.lll994f). We ca n estimate the 
ionization parameter of Kl. From the nuclear spectrum ([Wang et al.ll2010bl ). we obtain an 
ionizing luminosity Li 3 ^ eV _ W0ke y ~ 10 44 erg s -1 . Assuming a number density of 200 cm -3 
(Penston et al. 1990; Robinson et al. 1994) and a radial distance of 1.2 x 10 21 cm to the 
nucleus, we obtain an ionization parameter log U = —2. This is lower than, but in agreement 
with, the value of logU ~ —1.7 derived from the [OIII] /[Oil] ratio in Robinson et al. (1994), 
supporting the photoionization assumption. 

Since the X-ray blob and the [OIII] cloud are located at the same radial distance, their 
spatial offset suggests that the X-ray blob is in a higher ionization phase than the [OIII]- 
bright cloud. A likely explanation is that the X-ray blob has a lower density in the clump. 
If their densities are similar, the difference might arise from a higher ionizing flux received 
at the surface of the X-ray blob, perhaps due to less screening than in the direction towards 
the [OIII] cloud. 
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3.3. The Leaky "Torus" 

Our spectral fitting results show that a range of ionization states are needed to obtain 
good fits, in agreement with the multiple ionization parameters required in the grating 
spectra (e.g., Armentrout et al. 2007). 

A high ionization (logU ~ 1) component is found in the spectral fitting for all regions 
surrounding the nucleus, including the NW and SE sectors perpendicular to the ionization 
cone (Table 2). A low ionization (logC/ ~ —1) component is also present in the NW- 
SE direction that has lower U than found in the cone direction (NE-SW). This strongly 
indicates that there is no continuous absorbing "torus" blocking all the ionizing photons. 
This hypothesis is supported by the agreement between the observed H2 line flux in the 
"torus" direction and the value predicted from X-ray-irradiation by the active nucleus (Paper 
I), and by the measurement of HI absorption in the nuclear region, where high levels of 
clumping on the smallest scales is found (Mundell et al. 2003). We note that a "patchy 
torus " or clumpy wind arising fr om the accretion disk in which the clouds are optically thick 



(e.g., lElitzur fc Shlosmanl 120061 ). may allow leakage of nuclear continuum for such a high 
U component. However, it cannot reproduce the logU ~ — 1 component in the directions 
where the continuum is blocked. 

Extended emission along the "torus" direction is clearly present in the large scale X-ray 
images shown in Figures [2] (indicated by the arrows) and 15k. The base of this emission is 
at the location of the nucleus, thus may indicate ionization by the leaked photons in this 
direction. Moreover, we found an intriguing coincidence between radial regions of enhanced 
X-ray surface brightness and the P.A. of outflowing optical clouds identified by the HST 
study (Das et al. 2005). They pointed out that these faint clouds flowing close to the plane 
of the torus would need to be ionized by leaked nuclear radiation, which is attested by our 
X-ray data. The low U suggests a lower ionizing flux received by some clo uds, likely due to 



nucle ar continuum filtered by "warm absorbers" (e.g.. lKrongold et al.ll2007l ) or the cone walls 



(e.g., iKraemer et al.l 120081 ) . Evaluating the difference between the ionization parameters in 
the direction of torus and in the bicone, we estimate that approximately 1% of the total 
ionizing flux is seen by the NW-SE clouds. 



3.4. AGN Interaction with the Host Galaxy Disk 

Previous work on the NLR kinematics of NGC 4151 have established that the bi-cone 
geometry in NGC 4151 results in one side of the conical outflow exiting the galactic plane 
at a small angle, partly intersecting the host disk (e.g., Robinson et al. 1994; Crenshaw 



- 15 - 



et al. 2000; Das et al. 2005; Storchi-Bergmann et al. 2009). Thus part of the host ISM 
is likely illuminated by the AGN radiation. Our X-ray work finds that multiple phases of 
photoionized gas are present in the ENLR of NGC 4151. The fitting results for the "OVII" 
rich (red in Figure E]) and "NelX" rich (green in Figure E]) regions indicate that there is a 
range of ionizations, with a lower ionization phase associated with the clumpy clouds and 
the higher ionization with the diffuse outflow. There is a lowest ionization (log £7 ~ —2) 
component at larger radii in the SW and NE sectors. This component may come from 
high density material in the host disk plane, perhaps the nuclear spirals, illuminated by the 
nuclear radiation. 

We further note that interactions between the biconical outflow and the host galaxy 
ISM, similar to the radio jet-cloud interaction in the nuclear region, may contribute to the 
heating of any hot gas present in the host disk. Currently there is no evidence in the gas 
kinematics that supports any strong interaction (Kaiser et al. 2000). However, as we discuss 
below, there is some morphological evidence that, at the largest radii, the outflow runs into 
a clump in the inhomogeneous ISM and perhaps results in local shock heating of the gas. 

Figure [11] gives a large scale (1 arcminxl arcmin) view of the key ISM components in 
the circum-nuclear region of NGC 4151, illustrating the directions of the kinematic major 
axis of the host galaxy (green solid line; P.A.~22°, Pedlar et al. 1992; Mundell et al. 1999), 
the large scale "weak fat bar" (green dotted line; P.A.~130°; Mundell & Shone 1999), the 
ENLR bicone (cyan line; P.A.~65°, Evans et al. 1993), and the distribution of HI (Mundell 
et al. 1999) and CO (Dumas et al. 2010) gas, respectively. 

In Paper I we noted that the NE soft X-ray emission reaches as far as 12" (780 pc) from 
the nucleus, but that the surface brightness then decreases significantly by a factor of 5. 
Comparing the continuum-subtracted Ha emission (Knapen et al. 2004) with the extended 
soft (0.3-1 keV) X-ray emission (Figure IT7b). the ionized gas traced by bright Ha emission 
is mainly located in the r ~ 10" biconical region along the NE-SW direction centered on 
the nucleus, which also closely follows the soft X-ray emission. An arc-like feature in the Ha 
image (~ 15" x 2") is clearly present at the r = 10" location (Figure [T6l) at the edge of the 
X-ray emission, and the northern CO gas lane is also found here (Figure [TTb. and Dumas et 
al. 2010). 

Figure [TTb provides a multiwavelength, enlarged view of this region, showing the spatial 
relations between the sof t X-ray emission, the Ha, the CO gas lanes, and the HI distribution 
jMundell fc Shone! Il999h . The termination of the NE X-ray emission appears to be closely 
associated with the presence of the dense cold neutral material (traced by the CO and HI 
emission) in the galactic disk. No similar feature is seen to the SW of the nucleus. 
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The origin of the Ha arc (Figure H7h ) at the terminus of the NE X-ray emission is 
intriguing. It is spatially very close to the northern CO ga s lane and is one of the dust arcs 



driven by the bar potential (IMundell et al. 



identified in the optical V — I color map (Asif et al.lll998f). which a re gaseous compressions 



1999 



Dumas et al.ll2010f ). In this context, the Ha 



arc may be a more clumpy continuation of the dense CO arc and photoionized by the AGN 
radiation. The rate of ionizing photons from the AGN is approximately 10 53 s _1 (Kraemer et 



al. 2005; c.f. Li 3 Q eV _ W0keV ~ 10 erg s , Wang et al. 2010b). Assuming that the arc has a 
line-of-sight depth that is comparable to its projected length (~1 kpc), the arc's intercepted 
ionizing photon rate is 8 x 10 51 s _1 at its distance to the nucleus (970 pc), a covering factor 
of almost 10%. This photon budget should be treated as an upper limit as some fraction of 
the ionization photons will be absorbed or scattered on the way to the arc. If we assume 
that th e gas is in photo ionization equilibrium, for the case B effective recombination rate 
of Ha (jOster brockl Il989l ) . a rate of 8 x 10 49 s" 1 ionizing photons s" 1 is needed to produce 



the observed Ha luminosity of the arc (Lhq ~ 8 x 10 erg s ), which is only 1% of the 
maximal rate of available ionizing radiation in the direction of the arc (for a covering factor 
f c = 10%). The scenario where part of the CO gas lane is photoionized by the AGN and 
produces the Ha arc seems plausible with a modest covering factor of f c ~ 0.1%, implying 
that the arc has a filamentary structure with a thickness of a few pc. 

On the other hand, we consider the possibility that the arc could be a bow shock 
feature from interaction between the biconical outflow and dense gas in the host disk piled 
in the CO gas lanes by the large scale stellar bar. Assuming that the X-ray emission is 
due to shock heating as the outflow encounters the dense CO lane, and that the shock is 
strong and adiabatic, we estimate the shock velocity v s , adopting a postshock temperature 
T ps = (3 / 16) ([iv s/Icb) (Lehnert et al. 1999), where \i is the mass per particle (CO and HI) 
and ks the Boltzmann constant. For the kT ~ 0.3 keV X-ray temperature measured at 
the arc, the required v s is ~150 km s _1 , which is much higher than the local sound speed 
c s = 10 km s -1 in the ~ 10 4 K Ha-emitting gas. Although there is no direct measurement 
at the exact position of the Ha arc, the measured outflow velocity in [OIII] at r = 8" NE of 
the nucleus is fully consistent with this value of v s , e.g., cloud#28 which has v = 190 ± 32 
km s" 1 (Kaiser et al. 2000). 



3.5. Mass Outflow and Kinematic Power 

For two of the brighter, spatially resolved [OIII]-emitting clouds seen in the HST images 
(clouds #1 and #2 in Figure [TU see also Figure 4 in Kaiser et al. 2000), we were able to 
constrain the number density for the X-ray emitting gas with Cloudy modeling of their X- 
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ray spectra when an inner radius r = 100 pc to the nucleus and an ionizing photon rate 
Q = 10 53 are set. The best-fit density is logn# = 2.9 ± 0.5 cm -3 . We use this number 
density to estimate the mass loss rate and the kinematic power of the hot phase outflow, 
which are key quantities that gauge the impact of AGN outflow on the host galaxy ISM. 



Following Barbosa et al. (2009) and IStorchi-Bergmann et all (120101 ). we calculate the 



mass loss rate using M x = nnmnVrCgA, where v r is the radial velocity (taken as the approx- 
imate outflow velocity), m H the proton mass, C g the gas filling factor, and A the area in the 
cross section of the outflow. For a radial distance of r = 4 x 10 20 cm (130 pc) and the hollow 
cone geometry of Das et al. (2005), we obtain A = 7 rr 2 (sin 2 6 nnt — sin 2 6i n ) = 1.2 x 10 40 cm 2 . 



Adopting C g = 0.11 (IStorchi-Bergmann et al.l 120100 and v r ~ 750 km s 1 (Crenshaw et al 



2000; Storchi-Bergmann et al. 2010), we obtain a mass outflow rate M x = 2.1M yr _1 and 
a kinetic luminosity for the outflow L out fi ow = l/2M x v 2 ~ 1.7 x 10 41 ergs s _1 , or ~ 0.3% of 
the bolometric luminosity of the AGN (L& ; = 7.3 x 10 43 erg s _1 ; Kaspi et al. 2005). 

The M x derived from the hot phase outflow is comparable to the v alue derived from 



NIR studies of ionized gas in NGC 4151 (IStorchi-Bergmann et al.ll2010l ). which estimated 
mass outflow rate sal. 2 M Q yr _1 along each cone. The M x value is over 10 times higher than 
the outflow rate derived from UV and optical spectra (0.16 M Q yr _1 ; Crenshaw & Kraemer 
2007), but the latter is measured much closer to the nucleus (at 0.1 pc). The optical emitting 
gas likely has a much smaller filling factor than the X-ray emitting gas, which explains the 
difference in the measurements. If our study measured an outer part of the same outflow, it 
could have accelerated to a higher velocity or have loaded ISM between 0.1 pc and 100 pc, 
hence a higher M x . 

How does M x compare to the accretion rate of the SMBH in NGC 4151? The bolometric 
luminosity of NGC 4151, L bol = 7.3 x 10 43 erg s" 1 (Kaspi et al. 2005), is about 1.2% 
of its Eddington luminosity, which is L^m ~ 6 x 10 45 erg s _1 for the black hole mass 
M BH = 4.57±g;|? x 1O 7 M of NGC 4151 (Bentz et al. 2006). Adopting L bol = r]M accr c 2 
and rj = 0.1 for the efficiency of the accretion disk (Shakura & Sunyaev 1973), this implies 
that the central engine of NGC 4151 is accreting at M accr = O.O13M yr _1 . Since the 
measured outflow rate is 160 times the accretion rate necessary to feed the active nucleus, 
there must be either significant entrainment of the host galaxy's ISM while the nuclear 
outflow is expanding, or part of the gas fueling the nucleus is ejected before radiating, or a 
high fraction (over 99%) of the X-ray gas does not take part in the outflow and is illuminated 
by the AGN "in situ". However, the latter seems unlikely as the grating spectra strongly 
indicate that the emission lines have outflow velocities of 200-300 km s -1 (e.g., Armentrout 
et al. 2007). Significant mass loading of the outflow by the host ISM may be most plausible 
and common, since previous estimated NLR outflow rates in Seyferts (~0.1-10 M Q yr _1 ) 
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generally exceed tenfold the accretion rate (Veilleux et al. 2005). 

Nevertheless, even with such a high M x in NGC 4151, the kinematic power L out fi ow 
suggests that only 0.3% of the available accretion powe r is extracted to drive the outflow, 



similar to the findings of IStorchi-Bergmann et all (120101 ) and Holt et al. (2006, 2011). For a 



sample of Seyfert NLR outflows (NGC 4151 not included) studied in the [SIII] A9069 emission, 
Barbosa et al. (2009) reported a lower L out fi ow / L bol = 10 -5 — 10~ 4 . As pointed out in Mathur 
et al. (2009), this poses a problem to the majority of quasar feedback models which require a 
higher fraction of the accretion power of the black hole (~5% of L^z to unity) to be thermally 
coupled to the host ISM. Our L out fi ow / 'Lbol is i n rough agreement with a two-stage feedback 
model recently proposed by Hopkins & Elvis (2010), which requires only 0.5% of the radiated 
energy to drive the initial hot outflow to efficiently shut down star formation in the host. 



4. Conclusions 

In this paper we present spectral analysis and emission line images from deep Chandra 
observation of NGC 4151, aiming to resolve and characterise the X-ray emission in the NLR. 
The findings are summarized as follows: 

• The soft X-ray emission line images of NGC 4151 (OVII, OVIII, and NelX) are clearly 
extended and show remarkable morphological coincidence with the biconical NLR 
mapped by the [OIII] emission, which supports a common emission mechanism for 
the hot- and cool-phase of the NLR gas. 

• Extended emission in the X-ray image is detected along the NW-SE sectors, which is 
the direction of a putative torus. This may explain the faint rogue clouds identified 
in previous HST studies that require ionization in this direction, indicating leakage of 
nuclear ionization instead of full blocking by a continuous obscuring torus. 

• Spectral models involving smooth continua (a bremsstrahlung plus a power law) with 
emission lines provide good descriptions of the spectra. The emission lines cannot be 
uniquely identified with the present spectral resolution, but are consistent with the 
brighter lines seen in the Chandra HETGS and XMM-Newton RGS spectra below 
2 keV. The absorption corrected X-ray luminosity of the extended emission between 
r = 130 pc and r = 2 kpc is Lo.3-2kcV = 1.1 ± 0.2 x 10 40 erg s -1 . 

• Photoionization models successfully reproduce the soft X-ray emission, supporting the 
dominant role of nuclear photoionization. There are considerable variations in ioniza- 
tion states across the circum-nuclear region. A high ionization (log £7 ~ 1) component 
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is present in most regions. A low ionization (log U ~ —0.25) component is present along 
the bicone direction (NE-SW), whereas a lower ionization (log £7 ~ — 1) component is 
found in the NW-SE direction, which is consistent with filtered nuclear emission by 
warm absorbers instead of a continuous absorbing torus. The lowest ionization com- 
ponent (log U ~ —2) appears to be associated with the dense gas in the host plane. A 
thermal component may still be present at < 12% of the total soft emission, perhaps 
related to hot ISM in the galactic disk or shocks associated with the outflow. 

The measured ratios of [OIII]/soft X-ray flux are consistent with a constant ratio of 
~15 for the 1.5 kpc radius spanned by these measurements. This suggests a similar 
relative contributions from the low- and high-ionization gas phases at different radii. 
If the [OIII] and X-ray emission arise from a single photoionized medium, this further 
implies an outflow with a wind-like density profile (n e oc r~ 2 ). 

The estimated mass outflow rate in NGC 4151 is ~2M yr _1 at 130 pc and the kine- 
matic power of the ionized outflow is 1.7 x 10 41 erg s" 1 , 0.3% Ly. This value is signif- 
icantly lower than the expected efficiency in the ma jority of quasar f e edbac k models, 



but comparable to the two-stage model described in lHopkins fc Elvid (120101 ) . 



Placing all our findings in the context of previous studies, we obtain a comprehensive 
view of the NGC 4151 circum-nuclear region at various spatial scales, from the inner-most 
r ~ 50 pc (as illustrated in the schematic drawing in Figure H8k ) to as far as r ~ 2 kpc 
(Figure [18b). The key points are recapped here. 

Photoionization by the nucleus is important at all scales. In the nuclear region, except 
for a few clouds that show interaction with the radio jet (Wang et al. 2009, 2011b), most 
of the X-ray NLR clouds are consistent with being part of a photoionized biconical outflow, 
with ana r~ 2 density profile expected for a nuclear wind. 

In the ENLR of NGC 4151, photoionization is still the dominant ionization mechanism 
for the observed X-ray emission (e.g., Ogle et al. 2000; Kraemer et al. 2005), but a wide 
range of ionization states is present. This is likely related to the nature of the ENLR. The 
bi-cone geometry indicates that one side of the conical outflow exits the galactic plane at a 
small angle, partly intersecting the host disk (e.g., Crenshaw et al. 2000; Das et al. 2005; 
Storchi-Bergmann et al. 2010). Optical studies (e.g., Robinson et al. 1994; Kaiser et al. 
2000) have identified such a boundary between NLR and ENLR at r rj 6" (SW), and the 
ionized gas kinematics become consistent with the rotation velocity of the galactic disk at 
r > 6". Thus the origin of the ENLR is best explained as the inhomogeneous ISM in the 
galactic plane ionized by the AGN. Part of the high density nuclear spirals illuminated in the 
path of the bicone can produce the curvy shape of X-ray enhancement and the associated 
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low ionization. The presence of ionized emission perpendicular to the bicone indicates leaked 
nuclear emission along the putative torus direction, likely filtered nuclear emission by warm 
absorbers. 

In the large 3 kpc-scale cavity of the H I material, we find faint soft diffuse X-ray 
emission that provides evidence for AGN-host interaction in NGC 4151 (Wang et al. 2010a), 
originated from either hot gas heated by the nuclear outflow or photoionized gas from a recent 
nuclear outburst. 

As a concluding remark, our findings in NGC 4151 demonstrate that abundant in- 
formation can be extracted from X-ray spectral imaging studies when supplemented with 
multiwavelength data. Many valuable new radio, optical and IR observations and modeling 
of other nearby Seyfert galaxies are emerging (Crenshaw et al. 2010; Riffel et al. 2010; 
Fisher et al. 2011, Schnorr Miiller et al. 2011), and progress in the X-rays are being made to 
further our understanding of AGN feeding and feedback processes in these galaxies. 
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Fig. 1. — (a) Extraction regions of the large scale extended emission superposed on a gray- 
scale ACIS 0.3-2 keV band image. Excluded point sources are indicated by circles. The 
readout streaks running in the north-south direction is indicated, (b) The azimuthal surface 
brightness profile of the 0.3-2 keV emission, outlining sectors that contain bright extended 
emission. 
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Fig. 2. — (a) Adaptively smoothed ACIS image (0.3-2 keV) of the central 3 kpc-diameter 
region of NGC 4151. (b) Same as (a), but here the readout streaks are not removed, in 
order to avoid artificially increasing the significance of the NW-SE extension (indicated by 
the arrows). The contours are overlaid to help visualize the extent of the diffuse emission. 
The brightest nuclear region is masked. 
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Fig. 3. — (a) Comparison of the radial profiles of the bright SW cone (dots), the NW sector 
(asterisk), and the region in between (along P. A. ~ 270°; diamond), (b) Same as (a) but 
emphasizing the differences in the inner 5 arcsec region. 
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Fig. 4. — The residual spectrum of the large scale extended emission in NGC 4151 after the 
contribution from nuclear emission has been removed, demonstrating the presence of X-ray 
emission lines. Positions of the strong lines from the Ogle et al. (2000) HETG spectra are 
labelled, although they are blended at the spectral resolution of ACIS. The color bars outline 
the ranges used for narrow-band line images. 
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Fig. 5. — (a) Tri-color composite image of the central 3 kpc-diameter region of NGC 4151, 
where the soft (0.3-0.7 keV), medium (0.7-1 keV), and hard band (1-2 keV) adaptively 
smoothed images are shown in red, green, and blue, respectively. The inner r = 5" region 
is intentionally saturated and masked to show the faint features farther out. (b) Same as 
FigureE], but for the inner 1 kpc-radius region of NGC 4151. The inner r = 2" nuclear region 
is masked. 
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Fig. 5.— —Continued. 
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Fig. 6. — The narrow-band ACIS images of the three strongest emission features in the 
central 15" (~1 kpc-across) region of NGC 4151 (a) OVII; (b) OVIII Lya+OVII RRC; (c) 
NelX OVII triplet line emission. Overlaid are contours of the HST WFPC2 [OIII] A5007A 
image (Kaiser et al. 2000). The position of nucleus is indicated with a cross. "A" and 
"B" mark the X-ray structures displaced from optical clouds. "Kl" is a bright [OIII] cloud 
identified in Robinson et al. (1994). See § 13.21 (d) The radial profile comparison for the 
emission line fluxes with 1-a error bars, including the weaker emission lines NeX and Mg XI 
(see the legend in the top-left panel). 
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Fig. 6.— —Continued. 
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Fig. 7. — Spectrum of the extended X-ray emission with the best-fit phenomenological model 
consisting of a Bremsstrahlung component and gaussian lines, with lower panel showing the 
bin-to-bin contribution to the y 2 statistics. 
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Fig. 8. — (a) The same spectrum as shown in Figure 0, but with the best fit spectral model 
consisting of two photoionized components, (b)-(c) The photoionized models calculated by 
Cloudy for the low ionization parameter (logC/ = —0.25) and the high ionization parameter 
(log U = 0.8) component, respectively. 
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Fig. 9. — Comparison of the ACIS spectra across sectors and radially binned sub-regions 
(see Figure [Q, demonstrating the spectral differences, (a) the SW cone; (b) the NE cone; 
(c) the NW and SE sectors that are perpendicular to the bicone. 
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Fig. 10. — (a)-(d): Spectra of the X-ray emission in four regions radially separated along 
the SW sector (see Figure [Q . 




Fig. 11. — (a)-(c): Spectra of the X-ray emission in three regions radially along the NE 
sector. 
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Fig. 14. — HST [OIII] clouds labeled with numbers, for which the [OIII]/soft X-ray flux ratio 
were derived. 
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Fig. 15. — Radial dependence of the [OIII]/soft X-ray flux ratio. For comparison, the ratio 
for the Seyfert 2 galaxy NGC 7212 (Bianchi et al. 2006) is shown. 
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Fig. 16. — Continuum-subtracted Ha image of the central V x 1' of NGC 4151 obtained 
with the Jacobus Kapteyn Telescope (Knapen et al. 2004), which illustrates the directions 
of the kinematic major axis of the host galaxy (green solid line; P.A.~22°, Pedlar et al. 1992; 
Mundell et al. 1999), the large scale "weak fat bar" (green dotted line; P.A.~130°; Mundell 
& Shone 1999), and the ENLR bicone (cyan line; PA.~65°, Evans et al. 1993). The blue 
and red contours show the distribution of HI (Mundell et al. 1999) and CO (Dumas et al. 
2010) gas respectively. 
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Fig. 17. — (a) Superposition of a continuum-subtracted Ha image (contours; from Knapen 
et al. 2004) on the central 1' x V of the smoothed 0.3-1 keV ACIS image. The box region is 
enlarged in panel b. (b) A composite image showing the relative distribution of HI emission 
in red and Ha emission in green. Overlaid white contours are the soft X-ray emission, and 
blue contours the 12 CO emission. 
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Fig. 18. — A schematic drawing of the complex circum-nuclear environment of NGC 4151. 
(a) The inner few hundred pc radius region; (b) The 3 kpc-across region. The features are 
not drawn to exact spatial scale. The cyan wiggly lines represent possible leakage/scattering 
of ionizing photons, whereas the orange wedges denote the lower ionization sectors perpen- 
dicular to the bi-cone (yellow wedges). 
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Table 1. Blended Emission Lines Identified in the Large Scale Extended X-ray Emission 



ACIS 
Blend 


HETG 
Lines 


Observed Energy 3 - 
[keV] 


Rcstframc Energy 
[keV] 


Line Flux b 

[photon cur 2 s~ x ] 


ACIS/HETG Flux Ratio c 


OVII 


OVII f 
OVII i 
OVII r 


0.57 
0.5617 
0.5691 
0.5745 


0.5611 
0.5687 
0.5740 


5.5 X 10" 5 
3.1 x 10~ 4 
5.1 x 10~ 5 
6.8 x 10~ 5 


13% 


OVIII 


OVIII Lya 
NVII RRC 


0.66 
0.6541 
0.67 


0.6537 


1.7 x 10~ 5 
1.0 x 10~ 4 
2.6 x 10~ 5 


13% 


OVII 


OVII RRC 
OVIII Ly/3 


0.73 
0.74 
0.7749 


0.7747 


8.5 x 10~ 6 
5.2 x 10~ 5 
1.8 x 10~ 5 


12% 


OVIII 


OVIII Ly 7 
OVIII RRC 


0.82 
0.8175 
0.86 


0.8171 


4.4 x 10~ 6 
1.1 x 10~ 5 
1.1 x 10" 5 


2% 


NelX 


NelX f 
NelX i 
NelX r 


0.91 
0.9064 
0.9172 
0.9229 


0.9052 
0.9151 
0.9221 


1.3 x 10~ 5 
3.3 x 10~ 5 
9.9 x 10~ 6 
2.0 x 10~ 5 


21% 


NeX 


NeX Lya 
NelX ls3p - Is 2 


1.03 
1.0224 
1.0749 


1.0219 
1.0738 


4.4 x 10~ 6 
2.0 x 10~ 5 
5.9 x 10~ 6 


17% 


NeX 


NelX RRC 
NeX Ly/3 


1.20 
1.20 
1.2118 


1.2110 


1.9 x 10~ 6 

9.3 x 10- fi 

7.4 x 10~ 6 


12% 


MgXI 


MgXI f 
MgXI r 


1.30 
1.3542 
1.3525 


1.3313 
1.3524 


1.9 x 10~ 6 

1.2 x 10" 5 

9.3 x 10~ 6 


9% 


Sil 


Sil Ka 
MgXII Ly/3 


1.75 
1.7413 
1.7472 


1.7391 
1.7450 


1.6 x 10~ 6 
1.3 x 10" 5 
5.3 x 10~ 6 


9% 



Note. — The continuum is modeled with a brcmsstrahlung (best fit kT = 0.31 keV) absorbed by the Galactic column 
N H = 2.1 x 10 20 cm" 2 (Murphy et al. 1993). 

a Emission line energy measured in our ACIS spectrum and in Ogle ct al. (2000) HETGS spectrum. 

b ine flux from the present data and in Ogle et al. (2000) HETGS spectrum. 

c Ratio of the ACIS (present data) to HETGS line fluxes. 



Table 2. Spectral Fits to the Extended X-ray Emission 



Region 


Source 


logiVfl^.o.s 3 


logt/i 


logTV^i 


log{/2 


log N Hy2 


Nuclear PSF b 


X / d.o.f. 


Fo.5-2keV 




Counts [0.3-2 keV] 


[xlO 22 cm 2 ] 




[cm 2 ] 




r 9i 

[cm 2 ] 


Counts [0.3-2 keV] 




19 9 

erg s cm arcsec 


NE 1 


2432 


0.04±0.02 


0.86 


19.7 


-0.22 


19.7 


380 


63/65 


6.9 X 10" 15 


.... 2 


1814 


0.07±0.03 


1.05 


19.8 


0.06 


20.1 


200 


60/53 


1.3 x KT 15 


.... 3 


913 


0.03±0.01 


0.49 


19.9 


-1.87 


23.5 


140 


26/29 


1.7 x 10~ 16 


SW 1 


2096 


0.02 


0.83 


20.0 


-0.38 


19.4 


270 


60/57 


1.0 x 10" 14 


....2 


1667 


0.02 


1.06 


20.0 


-0.42 


19.8 


130 


58/45 


2.9 x 10" 15 


....3 


1560 


0.02 


1.00 


19.6 


-0.50 


19.6 


100 


54/46 


7.2 x 10- 16 


...A 


966 


0.02 


0.34 


19.0 


-1.28 


23.5 


80 


34/30 


4.5 x 10" 17 


NW 


1050 


0.22±0.02 


1.03 


19.2 


-0.98 


19.5 


400 


35/35 


3.1 x 10- 16 


SE 


1081 


0.16±0.04 


1.26 


20.3 


-1.16 


19.9 


530 


28/33 


2.2 x lO" 16 


OVII-rich 


2626 


0.02±0.02 


0.14 


20.3 


-1.22 


20.3 


200 


60/56 


4.3 x 10" 15 


NelX-rich 


2141 


0.07±0.02 


0.6 


19.6 


-0.21 


20.4 


210 


50/56 


4.0 x 10~ 15 



a The line-of-sight absorbing column density towards the photoioni zed emission. For t he SW regions, the fitted Njj i o a prefers values lower than the 
Galactic column towards NGC 4151, and was fixed at N H = 2 X 10 20 llMurphv et al.lll996h . 

b Expected counts from simulation of PSF scattered nuclear emission. 



Table 3. Measured X-ray and [OIII] Fluxes for the Clouds 



Cloud 


Distance 


Distance 


[OIII] flux 


0.5-2 kcV flux 


F[OIII]/ F[0.5-2keV] 


Number 


to Nuc. (") 


to Nuc. (pc) 


(10- 14 ergs s" 1 cm~ 2 ) 


(10 _15 ergs s _1 cm~ 2 ) 


1 


2.6 


169 


29.8±1.9 


25.7±5.2 


11.6 


2 


2.7 


175 


21.1±2.0 


15.2+1-° 


13.8 


3 


3.2 


208 


2.6±0.4 


9 A+°-3 
z '^-1.3 


10.7 


4 


4.8 


312 


7.9±0.9 


5 4+5-0 


14.5 


5 


3.5 


228 


1.2±0.4 


1 6 +0 - 2 


7.5 


6 


4.7 


306 


1.9±0.2 




4.1 


7 


6.1 


397 


10.2±6.0 




15.2 


8 


10.7 


696 


3.9±0.7 


4 2+°- 2 


9.6 


9 


17.9 


1164 


0.9±0.1 


2+°- 4 


43.4 


10 


17.7 


1151 


1.8±0.1 


nn+0.1 
u - M -0.8 


20.3 



Note. - 



The [OIII] flux is measured from F502N image. 



